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3 ■ ABSTRACT 
< 

We present timing analysis results for Rossi x-ray Timing Explorer (RXTE) 
observations of x-ray binary source 4U 1820—30 located in the globular cluster 

j> ■ NGC 6624. The light curves of observations made between October 1996 and 

September 1997 show that the maximum of the 685s binary period modulation 
folded by the linear ephemeris from previous observations has phase shift of 
—0.20 ± 0.06. Combined with historical results (1976-1997), the binary period 
derivative is measured to be P/P = (— 3.47 ± 1.48) x 10~ 8 yr _1 . The previous 
known (~176 d) long-term modulation is significant in the x-ray light curve 
obtained by analysis of the RXTE All-Sky Monitor (ASM) during the years 
1996-2000. The RXTE/ASM ephemeris is extended by analysis of all historical 

£j . data (Vela 5B and Ginga) to yield a period 171.033±0.326 days with no evidence 

for period change (\P/P\ < 2.20 x 10~ 4 yr _1 ). All reported x-ray burst activity 

£> . is confined to within ±23 d of the predicted minima. This stable long-term 

•i-H . 

^ , modulation is consistent with 4U 1820—30 being a hierarchical triple system 

with a ~1.1 d period companion. 

Subject headings: accretion, accretion disks — stars individual (4U 1820—30) - 
x-rays stars 



1. INTRODUCTION 

The low-mass x-ray binary (LMXB) 4U 1820—30 near the center of the globular cluster 
NGC6624 (Grindlay et al. 1984) has a binary orbital period of 685.0118 sec (Stella, Pried- 
horsky & White 1987; Smale, Mason & Mukai 1987; Morgan, Remillard & Garcia 1988; 
Sansom et al. 1989; Tan et al. 1991; van der Klis et al. 1993a,b) and was the first x-ray 
burster identified with a known x-ray source (Grindlay et al. 1976). Its short orbital period 
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and its x-ray burst activity imply it is a system with a 0.06-0.08 M Q helium white dwarf 
secondary star accreting mass onto a primary neutron star (Rappaport et al. 1987). A 
176 d period long-term variation between high- luminosity and low-luminosity by factor of 
~3 was observed by the Vela 5B spacecraft (Priedhorsky & Terrell 1984, hereafter PT84). 
Quasi-periodic oscillations (QPO) with various frequencies were also reported (Stella, White 
& Priedhorsky 1987; Hasinger & van der Klis 1989; Smale, Zhang & White 1997; Zhang et 
al. 1998; Wijnands, van der Klis & Rijkhorst 1999; Kaaret et al. 1999). 

Whereas the x-ray observations show an 11 minutes sinusoidal-like, small amplitude 
(~ 3% peak to peak) modulation, Anderson et al. (1997) discovered a large ~ 16% (peak 
to peak) modulation (period 687.6 ± 2.4 sec) in the UV band (wavelength in the 126-251 
nm range) from the Hubble Space Telescope (HST). The modulation may come from the 
variable thickness of the outer disk rim (Stella, Priedhorsky & Whitel987). The stability of 
the period (P/P = (-5.3 ± 1.1) x 10~ 8 yr -1 , van der Klis et al. 1993b) makes it certain 
that the 685 sec modulation is the orbital period. However, the negative period derivative 
is inconsistent with the lower limit (P/P > +8.8 x 10~ 8 per year) of the standard scenario 
proposed by Rappaport et al. (1987). 

The 4U 1820—30 luminosity variation was first discovered by Canizares and Neighbours 
(1975) and its possible ~176 d periodicity was first reported by PT84. The fact that x-ray 
bursts have only been seen in the low luminosity state (Clark et al. 1977, Stella, Kahn & 
Grindlay 1984) indicates that the long-term modulation is an intrinsic change in the accretion 
rate rather than an extrinsic absorption. However, significant phase shifts from the expected 
minima predicted by the ephemeris of the 176.4±1.3 d period proposed by PT84 have been 
observed from the EXOSAT (Haberl et al. 1987) and Ginga (Sansom et al. 1989). The ~176 
d periodicity, if stable, of the luminosity variation implies (Grindlay 1988) that 4U 1820—30 
may be a hierarchical triple in which a third companion star with a period of ~1.1 d orbits 
the 11 minute binary and thereby induces and an inner binary eccentricity precession (Mazeh 
and Shaham, 1979) with a period of ~176 d. 

In this paper, we describe our RXTE 1996-1997 PCA and 1996-2000 ASM observations 
of 4U 1820—30 (section 2) and report the timing analysis of the data (section 3), including 
the 11 minute binary periodicity, phase jitter, period stability, updated quadratic ephemeris, 
and search for the ~1.1 d period that could be associated with a third companion. Analysis 
of RXTE ASM data gives a new "176 d" modulation ephemeris which connects all the 
observations from 1969 to 2000 and exhibits no significant period derivative. In section 4, 
we discuss possible models for the observed light curve behavior of 4U 1820—30. 
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2. RXTE OBSERVATIONS 

The RXTE PCA/HEXTE pointed observations of 4U 1820-30 were made on 1996 
October 26, 28 and 30, and at least once per month between 1997 February and September. 
The observation time interval spanned about two 176 d luminosity cycles. Details of the 
RXTE PCA/HEXTE observations are listed in Table 2 in Bloser et al. 2000. The data used 
for the analysis are in PCA (PCU and 1) Standard-2 format with a time resolution of 16 
sec. We divide the data into 4 energy bands, 1.72-3.18 keV, 3.18-5.01 keV, 5.01-6.84 keV and 
6.84-19.84 keV. Since our timing analysis did not suggest any significant energy dependence 
of the patterns observed, we will present in this paper only the analysis results for band 4 
(6.84-19.84 keV) unless otherwise specified. A typical PCA Standard-2 light curve is shown 
in Figure 1. A complete analysis of the spectra and spectral variations of 4U 1820—30 for 
these RXTE observations is presented by Bloser et al. (2000). 

Data for the 4U 1820-30 RXTE/ ASM x-ray light curve (2-12 keV) analyzed in this 
paper were collected from 1996 January 11 to 2000 March 2 (Figure 2). The observation 
window (1576 days) spanned about 9 contiguous 176 d modulation cycles. In each modula- 
tion cycle, the count rate varies from ~5 to ~35 counts/sec. 

3. DATA ANALYSIS 

3.1. Binary Periodicity and Phase Analysis 

All RXTE PCA data were first corrected for the barycenter arrival times. In order to 
avoid the possible alias from the ~ 176 d period long-term modulation, we removed the DC 
term from the observed light curve. We carried out a x 2 analysis of the folded light curves 
(32 bins per period) to search for the best period near 685 sec for the entire 4U 1820—30 
RXTE PCA data. The maximum x 2 (deviation of folded light curve from constant flux) is 
at 685.014 sec as shown in Figure 3. Fitting the 685.014 sec peak with a Gaussian returns 
a best period of 685.0144 ± 0.0054 sec. The side bands are primarily artifacts of the alias 
period from the observation gaps. 

For comparison with historical results, we folded the light curve of each observation by 
the ephemeris from Tan et al. (1991). 

T N max = #71)2442803.63544 + (685.0118/86400) x N, (1) 
where N is the cycle count. The phase corresponds to the maximum of the sinusoidal fit of 
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each folded light curve. A typical folded light curve is shown as Figure 4. Figure 5 shows 
the phases of the RXTE 1996-1997 observations. The phase of maximum flux are scattered 
around -0.2 with ~ 0.061 (rms) phase jitter. 

To update the ephemeris, we appended the mean arrival times from the RXTE observa- 
tions to the historical results from SAS-3 (Morgen et al. 1988), Ariel V (Smale et al. 1987), 
Einstein (Morgan et al. 1988), Tenma (Sansom et al. 1989), EXOSAT 1 (Stella, Priedhorsky 
and White 1987), Ginga (Sansom et al. 1989, Tan et al. 1991 and van der Klis et al. 1993a) 
and ROSAT (van der Klis et al. 1993a,b). The phase (or arrival time offset) errors, however, 
need to be re-estimated. Van der Klis et al. (1993b) discovered the phase shifts about 0.038 
in the three ROSAT observations in 1991 and 1993, and the historical phase jitter ~0.050 
around the best fit linear ephemeris. For the RXTE data in this paper, we found the phase 
jitter is ~0.061 (=0.00048 d). Therefore, an additional 0.061 phase error was quadratically 
added to the historical data. For the RXTE data, we weighted-averaged the phases from 
each observation for the 1996 and 1997 datasets separately. We calculated the mean phases 
errors from quadratically adding the 0.061 phase jitter to the average phase errors from a 
sinusoidal fit of the folded light curves. The resulting mean phases are —0.20259 ± 0.0613 
and —0.20352 ± 0.0612 for the 1996 and 1997 datasets respectively. The mean arrival time 
for the average phases of two datasets were obtained from the expected flux maxima of the 
mid observation times of the observation windows, that is 



+ Pfold x {Nrnid + 'Pave) (2) 

where T = #JX>2442803.63544, P fold = (685.01180/86400) d, N mid is the cycle number 
closest to the mid time of the dataset and 4> ave is the mean phase for the 1996 and 1997 
observations. 

The period derivative can be obtained from a quadratic fit 

We applied linear (P = 0) and quadratic fits to the data. Both fits give acceptable results 
- x 2= 12.47 (d.o.f.=21) for linear fit and x 2= 6.96 (d.o.f.=20) for quadratic fit. However, the 
F-test (Bevington 1992) shows that F(z/!=Az/=l; z/ 2 =20)=15.83 for linear and quadratic fits. 
This implies that the quadratic fit is better than the linear fit at the ~99.99% confidence 



1 Arrival time errors of 0.0002 d were quadratically added, see van der Klis et al. (1991a) 
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level. Thus, the quadratic ephemeris is still required. 

From the quadratic fit of the data (Figure 6), we obtained P = (-7.54 ± 3.21) x 10~ 13 
sec/sec or P/P = (— 3.47 ± 1.48) x 1CT 8 yr" 1 , which is consistent with the value found by 
van der Kils et al. (1993b). The quadratic ephemeris 2 can be written as 



T N = #JX>(2442803.63564±2.2 x 10" 4 ) + (685.0119 ± 1.02 x 10~ 4 )/86400 x N 

+ (-2.99 ± 1.27) x 10~ 15 x N 2 (4) 

Tan et al. (1991) showed that the 685 sec modulation phases are well fitted by a period 
of ~8.5±0.2 yr sinusoidal curve from the 1976-1989 observation results. The period may 
be real or an artifact from the observation gap between 1981 and 1984 (Tan et al. 1991). 
To clarify the ambiguity, we used the constant-sinusoidal, linear-sinusoidal, and quadratic- 
sinusoidal models to fit the phases from all the 1976-1997 observations near the 8.5 yr period. 
The x 2 minimum fit results for the three different models are listed in Table 1. Although 
the linear-sinusoidal and the quadratic-sinusoidal models gave smaller reduced x 2 values 
than the quadratic model, the fitted amplitudes for both cases were only ~0.05 (modulation 
period ~6.5 yr), less than the 0.06 phase jitter. Therefore, there is no significant long-term 
phase periodic modulation of period ~6-8 years and the ~6.5 year "period" is highly likely 
to be an artifact from the phase jitter and the observation gaps between 1976 SAS-3 and 
1985 EXOSAT (-3500 d, 1.5 period), 1985 EXOSAT and 1989 Ginga (-1300 d, 0.5 period), 
and the 1989 Ginga and 1991 ROSAT (~1400 d, 0.5 period) observations. 

In explaining the discrepancy between the positive period derivative predicted by the 
standard scenario vs. the negative observed P, van der Klis et al. (1993b) demonstrated 
that the negative phase shift could be due to a long-term variation of the disk size. The 
proposal can be tested by looking for a dependence of orbital phase on M (van der Klis 
et al. 1993b), or L x . We compared the 11 minute modulation phases from the 1996-1997 
the RXTE observations and the simultaneous count rates from the RXTE All Sky Moni- 
tor (except February 9 1997 observation which has no All Sky Monitor data). The linear 
correlation coefficient is only 0.17, which implies that the uncorrelated probability is about 



2 The small offset due to leap seconds was ignored in our data analysis. The offset is about 16 seconds from 
zero phase epoch (JD2442803) to end of observations (September 9 1997=JD2450701). If we assume that 
the offset drifts linearly with time (first order approximation), this systematic effect will give only 1.6xl0~ 5 
sec offset in the second term of eq. (4). It is much smaller than the error (1.02xl0 -4 sec) from quadratic 
fit. The periodic systematic effect due to the difference between the heliocentric time and the barycentric 
time (~2.5 sec, mainly determined by the position of Jupiter) was also neglected. 
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70%. Therefore, no significant correlation between binary orbital phases and luminosities is 
observed. Further considerations about the phase shift are given in section 4. 

3.2. Possible Period Side Bands 

Grindlay (1986, 1988) suggested that the period of ~176 d long-term modulation of 
4U 1820—30 (PT84) may be due to a hierarchical triple companion star (captured by the 
compact binary in the high density cluster core) which modulates the eccentricity of the 
inner binary at a long-term period Pi ong = K P^ uter / P inner , where P inner and P ou ter refer 
to the binary period and the orbital period of the third companion and K is a constant 
of order unity which depends on mass ratios and relative inclinations (Mazah and Shaham 
1979). Under the triple model, with 176 d long-term modulation and 685 sec binary orbital 
period, the theoretical orbital period of the third star in the 4U 1820—30 system would be 
~1.1 d (for K~l, however, factors of 2 smaller or larger periods for the triple companion 
could be accommodated for differing inclination). If 4U 1820—30 is a triple system, the 
685 sec modulations would be affected by such a period and the beat side bands should 
appear near the peak of the power spectrum. We considered the binary motion around the 
center of mass of triple system. For the third companion star of mass ~0.5M Q (approximate 
maximum allowed by the optical counterpart) and ~1.1 d orbital period, the radius of the 
binary motion relative to the center of mass of triple system is only ~3.4 sim 3 light-second, 
where i% is the inclination angle of the orbit of triple companion. In other words, the observed 
~1.1 d period phase modulation amplitude is no more than 5xl0~ 3 . Although the phase 
variation form the binary motion may be too small to be observed, the third companion star 
could still affect the light curve in other ways. If, for example, the third companion star 
makes the 11 minute modulation amplitude change, ~1.1 d beat side bands 3 amplitudes may 
be detectable in the Fourier spectrum. 

To search possible side bands, we considered only the 1996 October 26 to 30 light 
curves because the observation gaps for the 1997 observations were too significant. To 
further minimize the observation windows (from observation gaps and Earth occultation), 
a one-dimensional CLEAN algorithm described by Roberts, Lehar & Dreher (1987) was 



3 The beat side band periods Pbeat = (1/ Pinner ± nj Pouter) 1 where n is positive integer and Pinner and 
Pouter are the binary (11 minutes) and third star period (<~1.1 d), respectively. The first harmonic (n=l) 
beat side band periods are thus 689.97 seconds (f=1.45xlCT 3 Hz) and 680.10 seconds (f=1.47xl0^ 3 Hz). 
The apparent sideband peak in the top plot of Figure 7 are tantalizing but probably due to the ~1 d spacing 
between successive observations. A longer continuous observation would be required to remove theses alias 
peaks. 
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applied to convert for the unequally spaced observations. We searched an arbitrary wide 
frequency range between 1.38 x 10~ 3 Hz (P=724.6 sec) and 1.55 x 10~ 3 Hz (P=645.2 sec) 
with amplitudes greater than 2a significance. The search results are shown as Figure 7. 
Only one primary peak is observed. The peak has an amplitude of 2.42 ± 0.26 cts/sec and 
a period of 685.120 sec. There are no other significant side bands beside this primary peak. 

3.3. ~176 Days Modulation 

The 1-day average binned light curve for the RXTE ASM data (see Figure 2) shows a 
clear modulation with relatively rapid rise and slow fall in count rate with a ~176 d period. 
The count rate varies from 5 cts/sec to 35 cts/sec (see Figure 2) throughout each cycle. 
"Inter dips" were also observed between lst-2nd, 4th-5th, 5th-6th and 6th-7th minima of 
the RXTE ASM light curve (minimum cycle count indices are marked on the bottom plot 
of Figure 2). 

PT84 reported the 4U 1820—30 long-term modulation period to be 176.4±1.3 d. How- 
ever, the ephemeris from PT84 (T min = JD2442014.5 + (176.4 ± 1.3) x N) does not match 
the EXOSAT 1985 August 19/20 4U 1820-30 low state observation (Haberl et al. 1985) 
where the expected minimum (by ephemeris from PT84) was offset by 48 days. A similar 
discrepancy was found in the data from the Ginga all sky monitor (delayed by ~50 d, San- 
som et al. 1989, S. Kitamoto, private communications). The light curve of the RXTE ASM 
data also shows a ~0.45 (80d) phase shift (see top plot of Figure 2). The inconsistency may 
be caused by an incorrect ephemeris (phase zero epoch, period or both) or a period drift. 

To obtain the best period to describe the 4U 1820-30 RXTE ASM light curve, we first 
applied a fast Fourier transformation (FFT). However, because the observation window is 
only a brief ~8.8 cycles (1576 days), the FFT frequency resolution, 5f = l/(1576d)= 0.233 
cycle per year, yields the period resolution near 176 d of 5P ~ P 2 x 5f= 19.8 d, which is too 
coarse for ~176 d period modulation. Therefore, the interpolated Fourier transformation 
(Middleditch, Geich and Kulkarni 1993) was applied to further determine the best period. 
The FFT is only able to show amplitudes of the "integer" frequencies (/„ = n/T where n is 
an integer and T is the total time of the data) whereas the interpolated Fourier may show 
the amplitude of "non-integer" frequencies (f r = r/T, where r is any real number). The 
non-integer amplitude A r can be estimated from the locally neighboring Fourier amplitudes 
as 
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[r]+m . , , ■ 

A r ~ V Aie ^r- ^n(n(r-l)) 

ir(r — I) 

l=[r]-m V ' 

where m is integer and [r] denotes the nearest integer of r. The uncertainty of peak frequency 



where 




3 d 2 P 

Po the peak power, and T is the length of the observation window (Middleditch, Geich and 
Kulkarni 1993). 

We chose m=2 and the resolution of r to be 0.1. The interpolated Fourier transformation 
spectrum is shown in the bottom plot of Figure 8. The peak amplitude was observed at 
f=2.130 cycle/yr with a value of 4.506 cts/sec. The frequency uncertainty from eq. (6) 
and eq. (7) is 0.0240 cycle/yr where the second derivative in eq. (7) is estimated by the 
quadratic fit around the peak of the power spectrum. Therefore, the best 4U 1820—30 long- 
term modulation period for the RXTE ASM data is 171.39 ± 1.93 d. Furthermore, the x 2 
period searching method gave the best period of 171.23 d (±7.36 d), close to the interpolated 
Fourier transformation result. 

The RXTE ASM light curve, derived by folding at a period of 171.39 d (see Figure 9), 
shows that the minimum closest to the mid of observation is at JD2450907.96 and ±0.07 
(rms) phase jitter (or ±12 d). The best linear ephemeris to describe the intensity minimum 
of the RXTE ASM light curve (hereafter local ephemeris) can be written as 



Tmin = ^2450907.96 ± 12.00 ± (171.39 ± 1.93) x N. (8) 

To obtain the best linear ephemeris for the intensity minimum in the RXTE ASM 
light curve and historical data, we assigned an uncertainty ±12 d (from 0.07 phase jit- 
ter obtained from the RXTE ASM light curve) to the minimum time reported by PT84 
(JD2442014.5). Combined with intensity minimum times observed by the Ginga all-sky 
monitor (JD2446822±26, S. Kitamoto, private communications) and RXTE ASM data, the 
linear fit yields the best historical ephemeris of 
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T min = JD2450909.90 ± 11.66 + (171.033 ± 0.326) x N. (9) 

The expected intensity minimum times from eq. (9) vs. the RXTE ASM light curve are 
shown as the bottom plot of Figure 2. An independent check on the ephemeris may be derived 
from timing of x-ray burst activity. Stella, Kahn and Grindlay (1984 and references therein) 
reported that no bursts are detected in the high state, which implies that the "burst phase" 
should be near phase zero. Table 2 lists the observation dates with bursts being detected 
and the phases of these days. No burst phase folded by eq. (9) exceeded the range ±0.13 
(±23 d). This indicates that the long-term modulation period of 4U 1820—30 is close to 
171.033 d and stable over ~30 years. 

To estimate the period derivative (or its upper limit), we suppose the minimum times 
(TV) obey a quadratic ephemeris (for small P) 

T N = T + P N + ^P PN 2 , (10) 

where N is cycle count index. Taking the phase zero epoch T to be that in eq. (8) 
(JD2450907.96), since the cycle numbers are small for RXTE intensity minima (from -4 
to 4), the quadratic term in eq. (10) can be neglected. Eq. (10) is reduced to a linear 
ephemeris equal to the local ephemeris (i.e. eq. (8)). The phases of historical data (Vela 5B 
and Ginga) folded by eq. (8) should be due to the period derivative 

where At is time difference between minimum time and To. Applying eq. (8) to historical 
data, we found no significant period derivative. The 2a (90%) confidence level upper limit 
for the change in 171 d period is \P\ < 1.03 x 10~ 4 days per day (= 0.038 days per year) or 
\P/P\ < 2.20 x 10~ 4 yr -1 . This stable long-term modulation is consistent with 4U 1820—30 
being a hierarchical triple system with a ~1.1 d period companion. 



4. DISCUSSION 

The observed negative period derivative from the 4U 1820—30 RXTE 1996-1997 ob- 
servations combined with all the historical data is consistent with the previous conclusions 
proposed by Tan et al. (1991) and van der Klis et al. (1993a,b). The decreasing 685 sec 
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period deviates from the positive P (P/P > 8.8 x 10~ 8 yr -1 ) predicted by the standard sce- 
nario (Rappaport et al. 1987). Tan et al. (1991) suggested that the discrepancy is probably 
caused by the acceleration of the binary star by the gravitational potential of the globular 
cluster. On the other hand, by analyzing the theoretical minimum of the period derivative, 
van der Klis et al. (1993a) found that the gravitational acceleration by the globular cluster 
is not enough to explain the observed results even if the line of sight is very close to the line 
connecting the binary with the center of the cluster at the projected separation 4±1 arcsec. 

However, King et al. (1993) measured the NGC6624 cluster center with the HST Faint 
Object Camera (FOC) and discovered that 4U 1820—30 is 0.66 arcsec from the cluster 
center. For an assumed distance of 6.4 kpc (Vacca et al. 1986, Haberl and Titarchuk 1995), 
this is equivalent to only 0.02 pc (projected) from the core. Using the model proposed 
by van der Klis et al. (1993a), the maximum gravitational acceleration along the line of 
sight could be a/c=2.5 x 10~ 15 sec -1 . Combining the period derivative derived in section 
3.1 and its value from the standard scenario, we determined the acceleration along the 
line of sight to be a/c=3.9 x 10~ 15 sec -1 , only ~50% larger than the maximum value. 
Therefore, given the uncertainties both in P and the cluster acceleration (i.e. center and 
mass model), gravitational acceleration by the globular cluster is still a possible explanation 
for the negative period derivative (also see King et al. 1993). 

Another potential explanation of the negative period derivative (or, negative phase shift) 
was proposed by van der Klis et al. (1993b). The 685 sec intensity modulation of the 4U 
1820—30, as for the dipping sources, is probably caused by the occultation by the accretion 
stream of the vertical structure at the edge of the accretion disk (Stella, Priedhorsky and 
White 1987; Morgan et al. 1988; Sansom et al. 1989; van der Klis et al. 1993b). The azimuth 
point of impact depends on the disk size. The bulge on the disk rim could shift by as much 
as ~ —120°, larger than the observed phase shifts of ~ —72° (—0.2 phase; the approximate 
value needed to account for the negative P). However, if we consider only the standard 
scenario, P/P ~ 8.8 x 10~ 8 yr~ x , the phase shift due to the period change is expected to be 
~ +0.75 from the 1976 SAS3 to the 1997 RXTE observations. This result implies that the 
total bulge phase shift is ~ — 0.95(— 350°). The bulge is unlikely to have such a large phase 
shift. Furthermore, the disk size should be highly correlated with the accretion rate and, 
of course, the luminosity L x . As discussed in section 3.1, no significant correlation between 
orbital phase and luminosity (171 d variation) is found in the RXTE 1996-1997 data. We 
hence conclude that the negative phase shift is unlikely to be caused by a variation of the 
disk size. 

In this paper, we derived the long-term modulation period and showed from the ephemeris 
(eq. (9)) that the 4U 1820—30 bursts are only observed in the low state. By re-analyzing the 
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historical data as well as tabulated burst activity time, we derived a period P=171.033±0.326 
d and a limit on \P\ < 1.03 x 10~ 4 days per day The high correlation between the burst ac- 
tivity and the luminosity suggests that the 171 d modulation is stable and indeed an intrinsic 
luminosity change rather than an extrinsic periodic obscuration. This luminosity modula- 
tion and (primarily) its long-term stability supports earlier suggestions for a hierarchical 
triple companion. The mass transfer rate is very sensitive to the Roche lobe radius, which is 
proportional to the inner binary separation. A hierarchical companion third star will induce 
an eccentricity variation in inner (11 minutes) binary with period Pi ong = ^ P^ uter / P inner 
(Mazeh and Shaham, 1979). When the minimum separation of the inner binary is small, the 
mass transfer rate and luminosity changes are enhanced. The triple model to 4U 1820—30 
system implies a ~1.1 d period third companion is responsible for the 171.033 d long-term 
intensity modulation. The triple companion star affects the orbital motion of the inner 
binary through beats of 685 sec binary period and ~1.1 d companion star orbital period. 
Our RXTE observations were not sensitive to this because of both data sampling (~1 d 
observation gaps) and the small amplitude expected. An additional test of the triple model 
could be conducted by a continuous, or optimally sampled long (> 3-10 d) observation of 
4U 1820—30 to measure the small (~3.4 sim 3 light-second, see section 3.2) phase shifts, or 
the possible modulation side bands (Figure 7) without 1 d sampling alias effects. 

The authors thank S. Kitamoto for invaluable assistant with the Ginga data, RXTE- 
GOF for help with RXTE data analysis, and the HEASARC for archival RXTE ASM data. 
This work was supported in part by NASA grant NAG5-3293 and NAG5-7393. 
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Fig. I.— Light curve of 4U 1820-30 (6.84-19.84 keV) observed by RXTE PCA on 1996 
October 28/29. 

Fig. 2.— The 4U 1820-30 light curve observed by RXTE All-Sky Monitor from 1996 January 
11 to 2000 March 2. The dashed lines are the expected minimum intensity times for the 
ephemeris from PT84 (top) and ephemeris of eq. (9) in section 3.3 (bottom). 

Fig. 3. — The result of x 2 folded period searching of RXTE PCA 1996-1997 observation. 
The Gaussian fit of the peak near 685.01 sec yields a best period of 685.0144 ± 0.0054 sec. 

Fig. 4. — The folded light curve of the observation on 1997 September 10 folded by eq. (1) 
and with the dc flux subtracted. The maximum of sinusoid fit is at about phase 0.8 (or -0.2). 
Inter-dip is observed at phase 0.5 to 0.6. 

Fig. 5.— The phases of 4U 1820-30 RXTE PCA 1996-1997 observation. The mean fluctu- 
ation of the phases is about 0.061. 

Fig. 6. — The phases of 4U 1820—30 11 minutes modulation from 1976-1997 observation 
folded by the ephemeris of eq (1). The dashed line represents the best quadratic fit result. 

Fig. 7. — The raw (top) and CLEAN (bottom) Fourier amplitude (i.e. absolute values of the 
Fourier transformation) spectra of RXTE 1996 October observations. The upper amplitude 
limit for the CLEAN spectrum is 0.52 cts/sec (2a). There is no clear side band beside the 
peak; therefore, the side bands in the raw spectrum are likely artifacts from ~1 d observation 
gaps. 

Fig. 8. — Raw FFT spectrum (top) and interpolated Fourier spectrum (bottom) of the 4U 
1820 RXTE ASM light curves. 



Fig. 9.— The 4U 1820-30 RXTE ASM light curve folded by the local ephemeris (eq. (8)). 
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Table 1. Phase modulation fit results 



Model 


Reduced x 2 


Period from 


Amplitude from 






sinusoidal fit (day) 


sinusoidal fit 


L a 


0.594 






Q b 


0.349 






c+s c 


0.782 


2978.3 ±140.7 


0.0965 


L+S d 


0.315 


2475.9 ±224.8 


0.0567 


Q+S c 


0.223 


2368.9 ±267.0 


0.0446 



a Linear model 
b Quadratic model 

c Constant ± sinusoidal model around P=8.5 year 
d Linear ± sinusoidal model around P=8.5 year 
c Quadratic ± sinusoidal model around P=8.5 year 
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Table 2. X-ray Burst phases 



Observation date 
(dd/mm/yy) 


Observation date 
(JD) 


Phase folded by 
PT ephemeris 


Phase folded by 
eq. (9) 


Reference 


18/05/75 


2442550.5 


0.039 


0.124 


1 


28/09/75 


2442683.5 


-0.207 


-0.098 


2 


15/03/76 a 


2442852.5 


-0.250 


-0.110 


3 


20/08/85 


2446297.5 


0.280 


0.032 


4 



a Mid observation date of 1976 March 11.5 to 19.5 



References. - - (1) Clark et al. 1976; (2) Grindlay et al. 1976; (3) Clark et al. 1977; 
(4) Haberl et al. 1987 
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